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Abstract. Until recently, wave-particle duality has been thought of as quantum principle 
without a counterpart in classical physics. This belief was challenged after (i) finding that 
average dynamics of a classical particle in strong inhomogeneous oscillating field resembles that 
of a quantum object and (ii) experimental discovery of "walkers" - macroscopic droplets that 
bounce on a vertically vibrating bath of the same fluid and can self-propell via interaction 
^ C"| with the surface waves they generate. This paper exposes a new family of objects that 

{"S j can display both particle and wave features all together while strictly obeying laws of the 

Newtonian mechanics. In contrast to the previously known duality examples in classical 
physics, oscillating field or constant inflow of energy are not required for their existence. These 
objects behave deterministically provided that all their degrees of freedom are known to an 
observer. If, however, some degrees of freedom are unknown, observer can describe such objects 
only probabilistically and they manifest weird features similar to that of quantum particles. 
We show new classical counterparts of such quantum phenomena as particle interference, 
tunneling, above-barrier reflection, trapping on top of a barrier, and spontaneous emission of 
radiation. In the light of these findings, we hypothesize that quantum mechanics may emerge 
as approximation from a more profound theory on a deeper level. 



1. Introduction 

Wave-particle duality is one of the fundamental principles of quantum mechanics which is directly 
linked to many of its mysteries. It postulates a dual nature of particles and radiation in 
microworld by enforcing the trade-off between the particle-like and the wave-like behavior of 
a quantum object, depending on experimental conditions. Since its original formulation as a 
hypothesis by Louis- Victor de Broglie in 1923 [1], the duality principle has been supported by 
numerous experiments involving photons [2], electrons [3-6], neutrons [7], atoms and dimers [8], 
small van der Waals clusters [9] and, more recently, Cqq fullerenes [10, 11]. Perhaps the most 
famous among these is the double-slit experiment, in which, let say, electrons are shot at a 
barrier with two slits, travel through the slits and are detected on a wall behind the barrier. 
The interference pattern formed by the electrons on the wall strongly suggests that they behave 
like waves in this setup. The interference pattern remains even when the electrons are shot 
separately, with interval between shots being large enough to ensure that only a single electron 
is present in the system at a time. If, however, an experimentator decides to determine through 
which slit each electron passes and sets up a detector close to one of the slits, mysteriously 
electrons stop creating an interference pattern, i.e., they start behaving like classical particles. 

Until quite recently, it was widely believed that wave-particle duality has peculiarly quantum 
nature, an attribute of behavior of matter at small length and energy scales that does not have 
a counterpart in classical physics. This belief remained unshaken by some known examples 



of 'classical tunneling' - crossings of potential barriers by objects of insufficient kinetic energy 
under condition that these objects are extended in space so that not all their parts are crossing 
the barrier region simultaneously (e.g., the Fosbury flop in high jump or a long train passing 
over the short hill) [12, 13]. The situation changed when in 2005 I. Dodin and N. Fisch showed 
that the average dynamics of a classical particle in a high-frequency oscillating field resembles 
motion of a quantum particle in a conservative field, with particle average displacement during 
the oscillation period playing the role of an effective de Broglie wavelength A. In a quasiclassical 
field, with a spatial scale large compared to A, the guiding-center motion is adiabatic. Otherwise, 
a particle exhibits quantized eigenstates in ponderomotive potential wells, tunnels through 
"classically forbidden" regions, and experiences stochastic reflection from attractive potentials 
[14-16]. In another development at approximately the same time, Y. Couder et. al discovered 
that macroscopic oil droplets can bounce on the vertically oscillating liquid bath and interact 
with the surface waves they generate to form stable self-propelled symbiotic structures called 
"walkers" [17, 18]. Walkers were shown experimentally to exhibit interference in double-slit 
experiment [19] along with other quantum-like phenomena such as unpredictable tunneling [20] 
and orbit quantization in a rotating frame [21]. 

The goal of this paper is to further extend the list of classical examples of wave-particle 
duality by presenting a rather general mechanism by which such examples can be produced 
in the framework of purely classical physics. Herein we present new classical counterparts of 
such quantum phenomena as particle interference, tunneling, above-barrier reflection, trapping 
on top of a potential barrier, and spontaneous emission of radiation. This list is likely to be 
continued increasing the number of links between classical and quantum theories. These links 
deserve serious attention and would be good topics for future research. Nevertheless, based 
on what is known so far including the findings presented in this paper, one cannot help but 
wonder whether the relationship between quantum and classical physics is much closer than is 
currently acknowledged. Especially intriguing is the question of whether quantum mechanics is 
an incomplete theory that emerges from the underlying classical sub-dynamics. This question 
has been the subject of long and on-going debates among physicists starting from the birth 
of quantum theory and has lead to many interesting theoretical developments including the 
stochastic mechanics introduced by E. Nelson [22, 23], statistical mechanics of matrix models of 
S. Adler [24], and emergent quantum mechanics based on nonequilibrium thermodynamics on a 
sub-quantum level pioneered by G. Grossing [25]. 

2. Classical mechanics of soft bodies 

Let us define a soft body (SB) as a collection of identical point particles (monads) of mass m 
coupling pairwisely with attractive central force. One can think of SB as a cloud consisting of 
interacting monads. Monads are confined to the cloud as a whole since parting the cloud into 
distant peaces would require a significant amount of energy. However, the cloud's overall shape 
can easily vary with time, thus distinguishing a soft body from a rigid one. The dynamics of 
monads inside SB is determined by coupling between them as well as by external forces that 
monads can also 'feel'. For simplicity sake, we assume that monad- monad interaction is given 
by an attractive harmonic potential, i.e., the force exerted upon a monad at location r-j by 
another monad at location Yj is proportional to a distance | r-j — r-j |. It follows that monads 
can penetrate through each other effortlessly (cf. asymptotic freedom) and, therefore, internal 
dynamics of monads is collisionless and can be treated analytically. Another assumption is 
concerned with how an external force acts on a monad. Let U(r, t) denote the external potential 
and N be a total number of monads in SB. We assume that a monad located at effectively 
'experiences' the external potential U(ri,t)/N. This calibration condition ensures that, when all 
monads happen to be co-located at the same point r, the total external potential of SB equals 



U(r,t). Thus, total potential energy of SB can be written as 
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where wo denotes the characteristic angular frequency of oscillations of SB. In what follows, 
we will be concerned with dynamics of the center of mass of SB moving in a smooth external 
potential in the non-relativistic limit. The main question is: How can we describe this motion 
effectively using as few degrees of freedom as possible? 

2.1. The simplest case: SB comprised from two monads in one dimension 
Consider the simplest possible instance of SB - two monads in one dimension. The total mass 
of SB reads M = 2m. Let x\ and x 2 denote the coordinates of two monads. Assume that 
the external potential is time- independent: U(x,t) = U(x). Also assume that, at any time 
instance, a linear size of SB is much smaller than the characteristic length at which U (x) changes 
substantially. It is convenient to introduce new variables as follows 

X = ^( Xl +x 2 ), £ = l(xi-x 2 ) (2) 

where A is a position of the center of mass (CM) and £ is an internal degree of freedom of 
SB. Expanding the potential U(x) in power series around CM, we obtain for the total potential 
energy Utot (up to the higher order terms) 

U to t(xi,x 2 ) = l - U( Xl ) + X - U(x 2 ) + l - mwg(xi - x 2 ) 2 U(X) + 1 Tj" (X)? + \ Muife 2 , (3) 

where prime denotes differentiation with respect to A. The total kinetic energy T tot of SB 
expressed via new variables reads 

T tot = X - MX 2 + l -Mi 2 , (4) 

where the overdot means differentiation with respect to time. From the action integral 

S = f\T tot - U tot )dt = f\\ MX 2 + \Mi 2 - U(X) - \ M(co 2 + ±u" (X))f }dt (5) 

we can then write down the Euler-Lagrange equations [26] as conditions for the action S to be 
stationary, by the Hamilton's principle, as follows 

MX = -U'(X)-^u'"(X)e (6) 

C + H + ^u"(xm = o (7) 

Notice rather remarkable relationship between the Eqs. (6) and (7). Equation (6) expresses 
the Newton's second law for the motion of CM, with the last term on the r.h.s. describing 
influence of the internal degree of freedom on this motion. Equation (7) can be easily identified 
as the time-dependent harmonic oscillator equation (TDHOE) for the internal degree of freedom 
£, with time dependence entering the picture via the term U [X(t)\. Thus, we see that CM 
parametrically affects oscillations of the internal degree of freedom, which in its turn, can directly 
influence the motion of CM (in the presence of position-dependent external potential). Also 
notice the mathematical equivalence of TDHOE to the time-independent Schrodinger equation, 
the root of interesting links between quantum and classical physics [27]. 



2.2. Conservation of energy 

When an external potential of SB is (explicitly) time- independent, the total energy is a conserved 
quantity. In this case, by combining (3) and (4), we obtain 

E = ^ MX 2 + U(X) + i Mi 2 + ^ M^ 2 + ^ ^"(^)£ 2 = constant (8) 



The first two terms in (8) give contribution that can be associated with the kinetic and 
potential energies of CM, the third and fourth terms add up to the 'internal energy', while the 
fifth term describes the coupling between (in this case) two degrees of freedom of SB and can 
result in energy exchange between them. Let us emphasize that, unless £ = or the external 
potential U is a constant or a linear function of X, the sum ^ MX 2 + U(X) does not conserve 
during motion. This is the main cause of quantum-like behavior of the CM of SB under various 
circumstances, as will be clear later. 



2.3. Time- dependent harmonic oscillator 
It is convenient to rewrite Eq.(7) as follows 

| + fj 2 (t)£ = (9) 

with Q 2 (t) = ujq + jj Jj"[X(t)]. The time-dependent parameter tt(t) is allowed to take not 
only real but pure imaginary values as well, so that f2 2 (i) may become negative for some time 
intervals without any contradiction [27]. Two linearly independent solutions of Eq.(9) can be 
written in the polar form 

£! = p(t)cos[8(t)} ; £ 2 = p(t)sin[9(t)} (10) 
and combined into one complex- valued function ip as follows 

<p(t) = Z 1 (t)+iUt)=P(t)e m (11) 



After substituting (11) into Eq.(10) and separately equating real and imaginary terms on 
both sides, we obtain two equations 

-|(P 2 0) = O (12) 
p at 

and 

p + Q 2 (t)p- p6 2 = (13) 



Equation (12) can be thought of as the angular momentum conservation law (for two- 
dimensional motion of a point particle of mass m in the plane (£1,^2)) and recast as 



^Mp 2 9 = J = constant, (14) 



where constant J has a dimension of action (or angular momentum) and plays an essential 
role in the theory to be developed below. Expressing 9 via p from (14) and substituting it into 
(13), we obtain the celebrated Ermakov-Pinney equation [28, 29] 

4 T 2 



Alternatively, one can assume that 9 > for all times and, from (14), find p 

and consequently p = —^J (9 ■ 9 — |# 2 )/# 5 / 2 . After substituting these results into (13), one 
finds that the phase function 9(t) must satisfy the equation 

n 2 (t) = 9 2 + ±{9;t}, (16) 

where {9;t} denotes the Schwarzian derivative [30] of 9(t) defined as {9;t} = 9/9 — |(#/#) 2 . 
Hence the (non-conserved) sum of kinetic and potential energies of the time-dependent harmonic 
oscillator, which coordinate is given by £i(t), can be written in the form 

ei(t) = \ Mil + \ MVi 2 ^ = Jco eff (t), (17) 

where the effective angular frequency ^ (t) is defined as 

u f At) = 9 + \cos 2 {9)-{9-9-9 2 )/9 z + \sin{29)-9/9\ (18) 

with 9(t) being the solution of Eq.(16) for a given Q 2 (t). 

The formulae (17) and (18) are general and therefore must be applicable in the particular 
case when fi 2 (i) = ojq = constant. In this case, the general solution of Eq.(16) is given by 

9{t) = arctan [k ■ tan (ujQt + ao)] (19) 

where k is an arbitrary positive number, and ao denotes the phase constant. The total energy 
of the oscillator is, of course, time-invariant in this case and reads 

ei = Jujq/k (20) 

In the light of (17) and (18), we can rewrite the energy conservation law for SB (8) as follows 

E = - MX 2 + U(X) + Ju Qff (t) = constant (21) 
2.4- Equation of motion for center of mass 

Assuming that loq 3> U"(X)/M for any X, we can employ the Wentzel-Kramers-Brillouin 
(WKB) method [31] to obtain the approximate solution of Eq.(7): 



J2J/M 

£(t) » . V ' ■ cos 

^l + U"[X(t)]/M 



(J y/<4 + U"[X(r)]/M dr + a) (22) 

where a is a phase constant, and we have used the condition that, in absence of the external 
potential (for free SB), £(i) = \/2J j '(Mujq) cos (u>ot + a) with a = a — woto- Notice the 
dependence of £(t) not only upon the value of X at moment t but also upon the whole path 
history X(t) at earlier times r < t via the integral in the argument of the cosine function in 
(22). Interestingly, this feature is reminiscent of "path memory" of walkers, which is believed 



to be responsible for quantization of their orbits in a rotating bath [21]. Substituting (22) into 
(6), we finally obtain the integro-differential equation of motion for a CM of SB: 



MX = -U\X) JU {X) = _ cQs 2 / /■* JuI + U"[X{t)]/M dr + a) (23) 

M^ujI + U"{X)/M Vito v - - y 

In case J = 0, i.e., when the internal degree of freedom of SB is not excited, Eq.(23) coincides 
with the second Newton's law for motion of a point particle of mass M in the external potential 
U(X). However, if J ^ 0, the behavior of CM can be significantly altered in comparison with 
that of a point particle of the same mass, presenting a great surprise to an observer who is 
unaware of that difference and treats SB as a point particle. It is worth emphasizing that 
Eq.(23) is valid even for a wider class of potentials that explicitly depend on time U(X, t), while 
the energy conservation law (expressed by (8) or (21)) does not hold in this case. 

To further simplify our analysis, henceforth we accept a crude approximation of (23) given 

MX = -U eff (X,t) (24) 
U et f (X, t) = U(X) + 3 '^W • cos 2 (wot + a) (25) 

where the phase constant a has been redefined to absorb — wq^o- 



3. Classical examples of wave-particle duality and analogues of quantum 
mechanical phenomena 

Here we present new examples of wave-particle duality in the realm of classical mechanics based 
on the theory developed in section 2. 



3.1. Classical analogue of tunneling 

Consider a potential barrier of a Gaussian shape U(X) = Aexp ^ y ~J ^ ^. Assume that a SB 
of mass M approaches this barrier from the left with CM moving at constant velocity V such 
that the corresponding kinetic energy KE = ^MV 2 (far from the barrier) is smaller than the 
barrier height A, as shown in Fig.l. Analysis of Eqs. (24) and (25) reveals that such SB may 
either reflect back or go through the barrier, depending on the value of the phase 9 = uj$t + a at 
moment when the CM is near the barrier top. Figure 2 shows the CM trajectories for different 
values of phase constant a chosen at random from the continuous uniform distribution on the 
interval [0,27r[. Thus, a plays the role of a hidden variable that completely determines the result 
of the interaction of SB with the potential barrier. For an observer who treats this SB as a point 
particle with the equivalent mass, its behavior seems unpredictable with apparent violation of 
the energy conservation law for some short time intervals and, therefore, can be classified as 
tunneling. 



3.2. Classical analogue of above-barrier reflection 

Consider now the case when KE > A as depicted in Fig.3. With appropriate choice of 
parameters, it is possible to obtain above-barrier reflection of a SB as shown in Fig.4. 



3.3. Classical analogue of particle trapping on top of a barrier and spontaneous emission 
Another example of paradoxical behavior of SB can be obtained if we select KE > A but 
KE — A <C A, see Fig.5. In this case, apart from the above-barrier reflection, we observe 
a new phenomenon: the CM can be confined for considerable interval oscillating around the 
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Figure 1. Setup for classical analogue of Figure 2. Trajectories for 100 different 
quantum tunneling. See text for details. values of phase constant a. 





Figure 3. Setup for classical analogue of 
above-barrier reflection. See text for details. 



Figure 4. Trajectories for 100 different 
values of phase constant a illustrating 
classical above-barrier reflection. 



top of the potential barrier and then be released in either direction with equal probability and 
approximately the same velocity as it had before the confinement. This behavior is similar to 
that of a quantum particle confined in a potential plateau region, as has been discussed in detail 
in [32]. Numerical simulations show that direction of the emission of a trapped SB is highly 
sensitive to the phase 6 of the internal degree of freedom, i.e., small change in the value of can 
lead to emission in the opposite direction. 

3.4- Interference in ID 

Interference phenomena are ubiquitous in Nature and are most easily understood using the 
theory of waves. For instance, colorful patterns produced by soap bubbles can be explained 
by constructive and destructive interference of light waves. It is much more demanding to 
explain thin- film interference using the representation of light as a collection of photons, requiring 
concepts of quantum theory, such as probability amplitude, to be invoked. Interestingly, Sir Isaac 
Newton - who considered light as being composed of tiny particles which he called 'corpuscles' - 
made some ingenious arguments in attempt to interpret the results of an experiment measuring 
the relationship between the thickness of a sheet of glass and partial reflection of light off it. 



>< 





Figure 5. Setup for classical analogues of 
paradoxical confinement of a particle on top of 
a potential barrier and spontaneous emission. 
See text for details. 



Figure 6. Trajectories for 100 different 
values of phase constant a illustrating 
classical confinement on top of a barrier and 
subsequent spontaneous emission. 



In the words of R. Feynman [33], "...to account for the fact that the thickness of the glass 
determines the amount of partial reflection, Newton proposed this idea: Light striking the first 
surface sets off a kind of wave or field that travels along with the light and predisposes it to 
reflect or not reflect off the second surface. He called this process "fits of easy reflection or easy 
transmission" that occur in cycles, depending on the thickness of the glass". 
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Figure 7. Scattering of a SB off a 'soft 
rectangular' potential barrier. See text for 
details. 



Under specific circumstances, a SB can exhibit behavior alike interference of waves in thin 
films, although it can be revealed only statistically as a build-up process of an interference-like 
pattern from single events and, besides, has a temporal character. Consider an ensemble of soft 
bodies of the same type prepared such that all members have fixed kinetic energy KE of their 
CM but different phase constants cc, drawn from the continuous uniform distribution on the 
interval [0, 2tt[. Each SB moves from the left toward the potential barrier of 'soft rectangular' 
shape U(X) = 4 (tanh(^t^ -tanh(^j^)) shown in Fig. 7. The width of the barrier is 
given by D = 2d. The kinetic energy of CM exceeds the barrier height: KE > A. We define 
the fraction of soft bodies in the ensemble trapped inside the barrier at some time instance 
t = T m as the trapping coefficient. The transmission coefficient is defined along the similar lines. 
Results of numerical simulations demonstrate that both trapping and transmission coefficients, 
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Figure 8. Trapping coefficient as a function 
of barrier width D estimated numerically 
after fixed time interval since the first collision 
of a SB with the barrier. 



Figure 9. Transmission coefficient as 
a function of barrier width D estimated 
numerically after fixed time interval since the 
first collision of a SB with the barrier. 



estimated at some fixed r m , depend quasi-periodically upon the width of barrier D, as shown in 
Figures 8 and 9. This oscillatory dependence fades away with increasing r m , as more and more 
soft bodies of the ensemble which are temporally confined inside the barrier start to leak out in 
either back or forward direction until none are left. 



4. Conclusions 

We have presented a new family of classical objects (soft bodies) obeying laws of the Newtonian 
mechanics that can exhibit particle and wave properties simultaneously. Such systems extend 
the list of known examples of the wave-particle duality in classical mechanics. They do not 
require an oscillating external field or constant inflow of energy to manifest their dual properties 
and behave deterministically provided that all degrees of freedom are known to an observer. If 
some degrees of freedom are unknown, soft bodies can be described only probabilistically and 
manifest weird features similar to that of quantum particles including interference, tunneling, 
above-barrier reflection, trapping on top of a barrier, and spontaneous emission of radiation. 
Further developments are underway. In the light of these findings, it is tempting to hypothesize 
that quantum mechanics may emerge as approximation from a more profound theory on a 
deeper level. If true, the road to achieving complete understanding of quantum theory may 
lie through an attempt to construct an average dynamics of classical objects like soft bodies, 
in which unobserved degrees of freedom have been effectively integrated out. This average 
dynamics will necessarily have a probabilistic nature. One can speculate that a concept of wave 
function may emerge as a mathematical tool to cope with lack of information about all degrees 
of freedom of a soft body, and the Schrodinger equation may even be deduced from the first 
principles. Such program is in line with the vision of A. Einstein who predicted: "Assuming the 
success of efforts to accomplish a complete physical description, the statistical quantum theory 
would, within the framework of future physics, take an approximately analogous position to the 
statistical mechanics within the framework of classical mechanics. I am rather firmly convinced 
that the development of theoretical physics will be of this type, but the path will be lengthy 
and difficult." [34]. The present paper advocates making new steps along this path. 
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